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PREFACE 


This  report  was  prepared  by  Hughes  Research  Laboratories,  Malibu, 
California,  under  contract  F30602-76-C-0022.  It  is  the  final  report  on  the 
contract  and  describes  the  analytic  and  experimental  algorithm  study  for 
thermal  blooming  compensation  using  a deformable  mirror  multidither 
COAT  system.  The  principal  investigator  and  principal  scientist  is 
Dr.  Richard  C.  Lind.  The  project  is  part  of  the  adaptive  optics  program  in 
the  Opto-Electronics  Department,  managed  by  Dr.  Viktor  Evtuhov,  at 
Hughes  Research  Laboratories. 
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SUMMARY  OF  PROGRAM  RESULTS 


The  objective  of  this  program  on  coherent  optical  adaptive  techniques 
(COAT)  is  to  investigate  algorithms  and  techniques  that  can  reduce  the  beam 
distortions  caused  by  thermal  blooming.  This  is  the  final  report  on  the 
contract  and  describes  the  analytic  and  experimental  work  performed  to 
determine  algorithms  that  can  compensate  for  thermal  blooming.  A descrip- 
tion of  the  previous  work  on  the  contract  can  be  found  in  Interim  Reports  1, 

2,  and  3 (Ref.  1,  2,  and  3). 

A major  result  from  the  work  conducted  on  this  contract  is  the  identi- 
fication of  certain  multidither  COAT  algorithms  employing  a deformable 
mirror  which  can  compensate  for  thermal  blooming.  The  specific  algo- 
rithms involved  and  the  degree  of  correction  achieved  depend  on  the  blooming 
scenario  studied.  Several  spatial  and  temporal  algorithms  were  studied 
both  analytically  and  experimentally  to  determine  their  effectiveness  for 
blooming  compensation.  In  particular,  under  the  analytic  task,  one- 
dimensional multidither  COAT,  function  maximization,  polynomial  phase- 
conjugate  algorithms,  and  initial  transmitter  irradiance  profile  algorithms 
have  been  numerically  investigated.  The  results  show  that  the  function 
maximization  algorithm  yields  the  greatest  correction  with  the  phase  con- 
jugate and  multidither  COAT  algorithms  yielding  somewhat  poorer  correc- 
tions. It  was  expected  that  the  multidither  algorithm  would  find  a solution 
intermediate  between  the  function-maximization  and  phase-conjugate  solu- 
tions, approaching  the  nonlinear  function-maximization  solution  as  the  servo 
convergence  time  was  slowed.  Instead,  as  the  servo  convergence  time 
(inverse  bandwidth)  varied  between  practical  limits  of  0.04  and  20  msec,  a 
linear  (nearly  phase-conjugate)  solution  was  retained.  We  also  found  that 
the  final  converged  multidither  solution  depended  on  the  initial  phase  front  of 
the  mirror.  Temporal  sequencing  of  various  polynomials  (in  this  case, 
Legendre  polynomials)  did  not  improve  the  correction.  However,  sequencing 
did  indicate  that  excellent  correction  could  be  achieved  using  only  one  par- 
ticular polynomial. 


Studies  performed  with  various  initial  transmitter  irradiance  profiles 
indicate  that  the  peak  target  irradiance  can  be  maximized  by  making  the 
aperture  distribution  as  uniform  as  possible.  This  conclusion  was  found  to 
hold  for  gaussian,  truncated  or  apodized  gaussian,  annular  or  circular 
beams,  and  for  special  hole-in-the-middle  laser  modes.  (See  Ref.  1 and  2.) 

Under  the  experimental  task,  zonal  multidither  and  modal  multidither 
employing  Zernike  polynomial  algorithms  were  investigated.  The  zonal 
multidither  achieved  varying  degrees  of  blooming  correction  over  a wide 
range  of  scenario  parameters.  The  modal  system  achieved  correction  for 
only  one  set  of  blooming  parameters.  We  speculate  that  the  modal  results 
suffered  because  only  lower-order  Zernikes  were  used  and  because  there 
was  cross  coupling  among  these  Zernikes.  This  problem  is  the  direct  result 
of  using  only  19  actuators  in  the  deformable  mirror,  which  resulted  in  the 
mirror  being  a poor  approximation  to  the  actual  analytic  Zernike  function. 
For  the  case  where  correction  was  achieved,  a sequencing  of  Zernikes  was 
performed.  It  was  found  that  no  improvement  in  correction  was  achieved  for 
the  various  sequences  performed  consistent  with  the  analytic  results.  In 
addition,  it  was  found  that  only  two  Zernikes,  focus  and  astigmatism,  were 
needed  to  give  the  correction,  again  consistent  with  the  analytic  results. 

Further  algorithms  employing  auxiliary  tracking  and  focus  servos 
were  also  studied  for  blooming  compensation.  Results  indicated  tracking 
worked  as  expected,  translating  the  wind-shifted  bloomed  beam  to  the  bore- 
sight,  while  focus  appeared  to  give  no  further  correction  over  what  the 
mirror  alone  would  give. 

As  a preliminary  to  the  use  of  the  deformable  mirror  in  the  blooming 
studies,  a characterization  of  the  mirror  was  performed.  A 37 -element 
all-beryllium  mirror  was  constructed;  however,  only  19  of  the  37  elements 
were  used  because  the  supplier  was  unable  to  deliver  all  of  the  mirror  actu- 
ators during  the  contract  period. 

Typical  performance  characteristics  of  the  mirror  were:  actuator 
surface  motion  of  0.2  pm  for  150  V drive;  resonant  frequencies  existing  in 
the  12  to  15,  18  to  20,  and  24  kHz  ranges;  an  influence  function  of 

1.71 

exp  (-5.396r  ' ) for  a faceplate  thickness  of  0.125  in.;  and  a mechanical 

coupling  coefficient  of  14%.  Convergence  times  of  1 . 5 msec  were  obtained 
and  tracking  limits  of  3 to  4 diffraction- limited  beam  diameters  were 
measured.  (See  Ref.  3.) 
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SECTION  1 


INTRODUCTION 

A.  PROGRAM  OBJECTIVES 

The  primary  objective  of  this  program  on  coherent  optical  adap- 
tive techniques  (COAT)  was  to  analyze  and  experimentally  demonstrate 
adaptive  multidither  correction  algorithms  which  can  reduce  beam  dis- 
tortions caused  by  thermal  blooming.  An  additional  program  goal  was 
to  investigate  the  use  of  fixed  transmitter  intensity  profiles  to  reduce 
thermal  blooming. 

B.  RESEARCH  PROGRAM  PLAN 

The  research  program  used  the  DARPA/RADC  experimental 
COAT  system  built  and  tested  on  contracts  F30602-73-C-0248  and 
F30602- 75-C-OOOl . Computer  simulation  codes  developed  on  these 
contracts  and  on  other  programs  (e.  g.  , NSWC  contract  N60921-74-C- 
0249)  were  used  for  the  analytical  portions  of  this  contract.  These 
codes  model  the  operation  of  several  types  of  COAT  servomechanisms 
and  model  the  time-dependent  propagation  of  optical  beams  in  an 
absorbing  and  turbulent  medium.  The  experimental  investigations 
required  constructing  a deformable  mirror  for  the  COAT  system.  To 
accomplish  the  contract  objectives,  a research  program  consisting  of 
four  major  tasks  was  developed.  Task  1,  an  analytic  task,  consisted 
of  computer  simulation  studies  of  blooming  compensation  algorithms. 
Task  2 was  to  design  and  construct  a 37-actuator  deformable  mirror. 

In  Task  3,  the  mirror  was  to  be  used  to  study  thermal  blooming  com- 
pensation with  zonal  multidither  control  and  with  modal  control  using 
Zernike-polynomial  modes.  Task  4,  an  amendment  to  the  contract, 
was  to  characterize  the  beryllium  deformable  mirror  that  was  built  and 
assess  the  applicability  of  the  design  to  high-power  laser  COAT 
systems. 
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REPORT  ORGANIZATION 


This  is  the  final  report  on  this  contract.  The  work  reported  previously 
in  three  interim  reports* will  not  be  repeated  here,  with  the  exception  of 
the  brief  "Summary  of  Program  Results"  presented  above.  The  interim 
reports  dealt  with  Task  2 (deformable  mirror  construction).  Task  4 (charac- 
terization of  the  mirror),  and  part  of  Task  1 (transmitter  profile  and  initial 
Zernike  blooming  studies).  The  primary  results  from  this  contract  which 
come  from  Task  1,  analytic  studies,  and  Task  3,  experimental  blooming 
compensation  studies,  are  reported  in  this  report.  Section  2 discusses  the 
analytic  work  performed  on  the  computer  to  study  various  algorithms  for 
blooming  compensation.  Section  3 describes  the  experimental  program 
conducted  for  blooming  compensation.  Section  4 discusses  the  Zernike 
polynomial  and  how  the  deformable  mirror  was  made  to  approximate  the 
polynomials.  Section  5 briefly  discusses  our  recommendations  for  future 
work. 
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SECTION  2 

ANALYTIC  BLOOMING  COMPENSATION 

The  goals  of  the  analytic  work  on  adaptive  algorithms  were  to  increase 
our  understanding  of  and  possibly  to  improve  the  dynamics  of  multidither 
response  to  thermal  blooming.  We  hoped  also  to  develop  a model  to  predict 
the  steady- state  behavior  of  a multidither  system. 

To  achieve  these  objectives,  two  areas  were  studied.  First,  we 
investigated  the  effect  of  varying  the  servo  bandwidth.  The  motivation  for 
this  study  resulted  from  the  fact  that  the  best  correction  for  thermal  bloom- 
ing was  invariably  predicted  by  a nonlinear  scheme  which  resembled 
infinitely  slow  multidither.  The  objective  was  to  move  toward  this  non- 
linear correction  by  slowing  the  servo.  Second,  we  analyzed  temporal 
sequencing  of  the  polynomial  phase-correction  channels  to  improve  the 
blooming  compensation.  Polynomial  multidither  systems  permit  channel 
signal  mixing  when  the  phase  distortions  are  large.  Thus,  simultaneous 
channel  convergence  of  a system  with  initially  large  errors  could  easily 
find  a nonoptimum  correction  path.  Sequential  channel  convergence  might 
yield  path  changes  significant  enough  to  find  a different  steady- state  solution. 

This  section  begins  with  an  examination  of  the  mathematical  models 
used  in  this  study,  including  the  choice  of  a one -dimensional  (1-D)  propaga- 
tion code,  the  linearities  and  nonlinearities  of  polynomial  phase  correction, 
and  a discussion  of  computer  codes  which  model  the  various  optimization 
schemes.  This  is  followed  by  the  results  of  the  bandwidth,  sequencing,  and 
predictive-modeling  studies  in  one  dimension  and  by  the  limited  corrobora- 
tive results  in  two  dimensions. 

A.  MODELS 

1 . One-Dimensional  versus  Two-Dimensional  Propagation  Codes 

Our  studies  used  the  first  five  Legendre  polynomials  to  correct  a 
simulated  one-dimensional  blooming  medium.  (Our  nonlinear  propagation 
codes  proceed  in  ten  steps  from  transmitter  to  focus.  A one- 
dimensional code  uses  line  rather  than  plane  apertures  at  each  propagation 
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step.)  We  chose  a 1-D  over  a two-dimensional  (2-D)  code  for  the  bulk  of 
the  work  to  allow  more  extensive  runs:  using  a 2-D  simulation  to  a steady 
state  is  prohibitively  expensive.  Using  a 1-D  approach  has  the  disadvantage 
that  a quantitative  correspondence  between  the  1-D  and  2-D  simulations  is 
lacking.  From  the  results  obtained,  blooming  compensation  is  better  in 
1-D  than  in  2-D.  In  addition,  the  Legendre  polynomials  have  no  exact 
analogs  in  Zernike  polynomials.  However,  we  believe  that  the  qualitative 
results  remain  valid. 

2.  Polynomials 


The  first  five  Legendre  polynomials  and  some  associated  Zernike 
functions  are  listed  below.  Note  the  difficulty  in  making  simple  analogs. 


Z j i = pcosG  = x 

Z20  ’ 2"2  ‘ 1 

Z22  = pZ  cos  26  = X2  - y2 

Z^j  = (3  - 26)  cos  6 

3 

Zjj  = P cos  36 
Z4Q  = 6p4  - bp“  + 1 
Z 4£  = (4p4  - 3 pZ ) cos  26 

4 

Z44  = P cos  46 
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P.  = -i(63x5  - 7 Ox3  + 15x) 

D o 


Z51  = ( 1 Op5  - 12p3  + 3p) cose 
Z53  = (5p5  - 4 p3)  cos  3G 


= p3  cos  50 


3.  Linear  and  Nonlinear  Polynomial  Response 

Any  system  which  deals  with  thermal  blooming  will  be  nonlinear. 
However,  in  this  report  we  will  consider  a linear  system  to  be  one  in 
which  phase-correction  measurements  are  made  instantaneously:  the 

medium  has  no  time  to  respond  to  trial  perturbations,  and  the  system 
is  optimized  on  a linearly  propagated  beam.  A nonlinear  response 
system  is  one  in  which  optimization  measurements  on  trial  perturbations 
are  made  only  under  steady-state  conditions. 

The  two  may  arrive  at  different  solutions.  Consider  a mirror 
surface  in  the  x plane  focusing  down  to  a point  xq.  The  focal  spot 
intensity  is 


I = U(xo)  U (xQ) 


(2) 


To  maximize  I with  respect  to  a particular  channel  coefficient  p,  requires 


5T- *°*J-  Re  ("V'TST 


(4) 


The  last  derivative  determines  linearity  or  nonlinearity.  Since  a linear 
system  does  not  allow  the  medium  to  respond  to  the  perturbation  in  p^,  it 
yields 


au(xo) 


gp — = I U(x)  (-jkPk(x))  (G(x,xQ)dx 


(5) 


A nonlinear  system  waits  for  the  Green's  function  to  respond,  yielding 


3U*(xo) 

~WZ~ 


1 

aG*(x,xn) 

/ U*(x) 

(+jkPk(x))  G.(x2x0)  t — - 0 

k 

dx  . (6) 


Now,  to  find  an  expression  for  the  correction  signal,  we  make  the 
following  simplifications: 

• Start  with  a uniform  beam:  u(x)  = 1 . 

• Ajsume  the  Green's  function  has  only  a phase 
distribution,  which  depends  on  the  phase-correction 
set  p = (pr  P2,  p3  . . . . ): 


G(x,xq)  = exp^-jk  ^ an(p)  Pn(x)). 

n 


(7) 
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This  yields 


jjp  = 0 = 2 • Rc  [ j U(x)  G(x)  dx  |U*(x)  G*(x)(jk) 


Eaa  (P)  I 

Pn<*>  ~W~  dX  * (8) 


Further,  assume  that  all  net  phase  perturbations  are  small,  i.e., 


k(Pn  + an(P,)  < 1 


This  allows  linear  expansion  of  the  exponentials.  Using  the  orthogon- 
ality of  the  polynomials,  we  find 


*1  J(Pk  + ak(P))  {Pn  + an((3))  8an(P)) 

apk*'  I gk  2-  i;  -apri 


= o , (9) 


where  depends  on  the  type  of  polynomial  involved. 


Legendre:  g 


Un  + i; 
2 


Zernike:  gn  = (n  + 1),  m = 0 


gn  = 2(n  + 1),  m i 0 . 


For  a linear  system  the  partial  derivatives  of  Eq . (9)  will  be  zero 
auch  that 
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Pk  = -ak(P) 


This  is  the  phase-conjugate  result.  We  will  examine  the  nonlinear  result  in 
more  detail  later  in  this  section  to  clarify  some  of  the  computer  results. 


4.  Function  Maximization  and  Polynomial  Phase-Conjugate  Models 

Function  maximization  (FM)  and  polynomial  phase  conjugate  (PPC) 
are  two  idealized  systems  which  give  corrections  that  lie  at  the  extremes 
of  nonlinearity  and  linearity.  Function  maximization  finds  the  nonlinear 
solution.  With  each  trial  mirror  deformation,  a new  steady-state  medium 
is  determined  before  glint  intensity  is  calculated.  One  polynomial  at  a 
time  is  optimized,  then  held  at  its  optimum  while  the  next  is  varied.  After 
one  maximization  cycle,  all  coefficients  are  varied  simultaneously  a few 
times,  then  the  sample-and-hold  sequencing  starts  again. 

Phase  conjugate  is  the  simplest  example  of  a linear  correction 
system.  Polynomial  phase  conjugate  uses  the  return  wave  phase  informa- 
tion of  an  ordinary  phase- conjugate  system  and  decomposes  it  into  a sum 
of  orthogonal  polymials.  The  new  outgoing  wavefront  is  then  assembled 
— from  only  those  polynomials  used  in  the  multidither  system  being  modeled. 


5.  Multi  dither  Model 


A slow  system  in  which  the  medium  can  come  to  steady  state  after  each 
dither  "step"  — (t^  » t^)  — should  approach  a nonlinear  solution.  A fast 
system  (Tg  <<  t^)  should  approach  a linear  solution. 

b.  Constraints  on  the  Servo  Parameters 

We  have  scaled  the  entire  servo  structure,  including  the  dither 
frequencies,  to  the  bandwidth  f . This  allows  the  interplay  between  Tg 
and  to  be  isolated  as  the  only  temporal  variable.  The  servo  scaling 
parameters  are  shown  in  Figure  1.  Most  values  have  been  empirically 
determined. 


c.  Dither  Period 

A constraint  exists  which  links  to  the  servo  bandwidth: 


D > 10  • 2irf 
o 8 

AD  > 4 • 2-rrf 

s 

< IT 

Td  “ D + NAD 
o 


empirical  experimental 
constraints 


sampling  theorem. 


f = servo  bandwidth 
s 

Dq  = lowest  dither  frequency 
AD  = distance  between  adjacent  dithers 
N = number  channels 


(12) 
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d.  Amplifier  Gain 

Each  type  of  system  (modal  or  zonal)  requires  a different 

calculation  for  the  amplifier  gain,  G,  so  that  it  remains  consistent 

with  the  desired  servo  bandwidth. 

We  need  first  to  look  at  the  open-loop  gain  g,  of  which  G is  a 

factor.  The  transfer  function  of  the  servo  is  dominated  by  a low-pass 

filter  with  a corner  at  f : 

o 

2(0  = (T-Mi/f0)  ' <l3) 

Bandwidth  f is  defined  as  the  point  of  unity  gain: 

|Z(fs)|  = 1 

or 


I 


Open-loop  gain  is  given  by 


g 


(14) 


(15) 


where  I /lQ  is  the  signal  at  frequency  w after  the  automatic  gain  control 
(AGC)  and  the  factor  2 comes  in  after  synchronous  detection.  In  an  N- 
channel  zonal  system,  with  a dither  amplitude  4*  ^ on  each  channel, 

I~  = 4 Jo  W,d)  Jlw,d)  (i  ““z)  * "IT  • (16) 
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which  gives 


(17) 


In  a polynomial  system,  the  signal  gain  varies  with  the  order  of  the 
polynomial. 


I 

w 

I 

o 


(18) 


where  was  defined  in  Eq.  (10)  for  both  Legendre  and  Zernike  poly- 
nomials. Uniform  channel  gain  is  ensured  by  defining  a variable  dither 
amplitude  in  each  channel. 


*d(n)  = Vn’ 

This  yields 


(19) 


e . Medium  Response  to  Dithers 

The  standard  Hughes  multidither  code  assumes  a servo  fast 
enough  that  the  medium  cannot  respond  to  the  dithers  (i.e. , a linear  cor- 
rection system).  It  does  this  by  heating  the  medium  only  at  discrete 
medium  update  times,  with  a beam  phase  averaged  over  the  dither  per- 
turbations. In  the  intervals  between  update  times,  a Green's  function  is 
used  to  linearly  propagate  the  dither  beam.  To  see  any  nonlinearity  in  the 
correction,  it  is  necessary  to  propagate  the  beam  nonlinearity  at  each 
dither  sample  time,  with  an  instantaneous  rather  than  an  averaged 
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phase  profile.  We  implemented  this  as  a code  option:  multidither  with 
medium  response  (MDMR).  This  option  led  to  expensive  code  runtimes; 
consequently,  only  a few  runs  were  made. 


B.  SERVO  BANDWIDTH  VARIATION 
1 . Results 

Multidither  code  results  were  expected  to  give  a nonlinear  solution 
with  a slow  servo,  and  a linear  solution  with  a fast  servo.  However,  as 
servo  convergence  time  (inverse  bandwidth)  varied  between  practical  limits 
of  0.04  and  20  msec,  we  retained  a linear  solution.  From  the  same  initial 
state,  the  system  in  all  cases  converged  to  the  same  position,  or  oscillated 
between  the  same  limits.  In  addition,  the  1-D  correction  was  greater  than 
that  found  in  2-D  thermal  blooming  compensation. 

Of  the  many  cases  studied,  two  are  of  primary  interest  (see  Table  1 
and  Figures  2 and  3),  both  with  uncorrected  Strehl  ratios  around  0.5. 


Table  1.  Propagation  Parameters 


Case  (I) 

Case  (II) 

a = 2 x 10'^/cm 

a = 2 x 10'^/cm 

P = 20  kW/cm 

P = 20  kW/cm 

v = 500  cm/ sec 

v = 127  5 cm/sec 

f = 2 km 

f = 2 km 

ft  = 0.02  rad/sec 

ft  = 0.005  rad/sec 

APO  = 22  J/m^ 

APO  =17.3  J/m3 

p =70  cm 
o 

p = 3b  cm 

ro 

P,/e  = 53.5 

Pl/e  = 26.75 
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Case  (I)  converged  steadily  to  a value  of  0.  9 for  all  servo  convergence 
times.  Even  at  Tg  = 20  msec,  and  allowing  the  medium  to  respond  to  the 
dithers,  no  difference  in  Strehl  ratio  was  observed. 

Case  (II)  had  less  slewing  than  case  (I).  In  the  time  allowed  for  the 
runs,  peak  intensity  oscillated  from  0.7  to  over  1.  The  oscillation  period 
was  approximately  15  msec,  the  time  it  took  the  medium  to  reach  a steady 
state.  In  this  case,  the  addition  of  medium  response  appeared  to  create  an 
instability;  the  good  correction  suddenly  deteriorated  at  58  msec  and  again 
at  85  msec  (not  shown). 


Discussion 


a.  Multidither  versus  Function  Maximization 

How  slow  would  the  servo  system  have  to  be  to  approach  the 
function  maximization  results?  The  requirements  for  expected  nonlinear 
correction  are  given  in  2.  A.  5 above. 


T < T < S£-i-= 

m ~ d fg(4(N-l)  + 10)  * 


Tm  is  now  found  to  be  approximately  15  msec,  which,  for  a 5-channel 


system,  gives 


f < 1.  3 Hz 
s 


t > 0.  8 sec 

s 


This  is  too  slow  for  any  practical  dither  system. 


b.  Multidither  versus  Polynomial  Phase  Conjugate 

Referring  to  Table  2,  we  see  that  multidither  consistently  found 
a somewhat  poorer  solution  than  did  PPC  when  starting  from  a flat  mirror. 

A better  maximum  was  found  when  starting  near  the  phase  conjugate  solution. 
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Table  2.  Multidither-PPC  Results 


System  (Case(I)) 

Target 

Intensity 

Polynomial  Solution 

Px(x) 

P2(x) 

P3(x) 

P4(x) 

p5(x) 

Phase  Conjugate 

0.  93 

-7.  5 

1 

© 

o 

2.9 

1 

o 

LP 

-1.  86 

Multidither 

(No  initial  correction) 

0.  9 

-9.8 

0.  1 

0.  02 

-0.  1 

-3.  54 

Multidither 
(Starting  at  phase 
conjugate  solution) 

0.  92 

-7.  8 

0.  08 

3.4 

-0.  3 

-2.  15 

The  better  solution  using  multidither  can  be  understood  with  reference 
to  Eqs.  (8)  and  (9).  Starting  near  the  PPC  coefficients  results  in  smaller 
initial  net  phase  distortions.  Thus,  the  exponential  linearization  is  more 
valid,  complex  interactions  between  propagating  polynomials  are  of  less 
importance,  and  the  system  is  less  likely  to  head  for  a secondary  maxi- 
mum. In  our  1-D  multidither,  however,  it  appears  that  the  secondary 
maximum  is  barely  lower  than  the  primary. 

The  PPC  code  gives  instantaneous  correction  to  the  medium.  The 
infinite  bandwidth  allows  some  initial  instability  in  the  correction  for 
case  II,  which  does  not  occur  for  multidither.  However,  the  two  settle 
down  subsequently  to  almost  identical  temporal  behavior. 

As  an  interesting  sidelight,  we  might  speculate  as  to  why  the  1-D 
correction  is  so  good.  Figure  4 sketches  the  focal-plane  patterns  of  the 
uncorrected  and  nonlinearly  corrected  beams  of  case  II.  For  comparison, 
we  have  included  the  focal  pattern  of  an  unbloomed  beam.  For  2-D  beam 
cross  sections  in  general,  it  is  unusual  to  see  a bloomed  beam  which  cor- 
rects to  a Strehl  ratio  greater  than  1,  but  not  impossible.  Wind  probably 
causes  the  worst  of  the  blooming  to  occur  some  distance  back  from  the 
focus.  If  the  beam  reaches  the  end  of  the  heated  medium  with  a somewhat 
wider  diameter  than  the  unbloomed  beam,  it  may  focus  down  through  the 
wind-cooled  section  with  a smaller  f#  to  a smaller  focal  spot.  However, 
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since  this  description  could  be  just  as  valid  in  2-D  as  in  1-D,  we  are  left 
with  the  aame  question:  Why  is  1-D  correction  so  much  better?  Reexamina- 
tion of  Lq.  (9),  which  delineates  the  linear  and  nonlinear  portions  of  the 
correction,  yields 


Two  points  should  be  noted: 

• This  equation  is  derived  from  a treatment  which  assumed 
only  phase  distortions,  and  therefore  has  limited  applicabil- 
ity to  a real  blooming  situation. 

• The  N terms  must  include  all  polynomial  coefficients  which 
make  a contribution  to  the  phase  description  of  the  medium, 
even  if  many  fewer  correction  channels  are  actually  used. 

It  is  easy  to  show  that  good  steady-state  correction  is  possible  if 
the  sum  of  the  nonlinear  terms  of  Eq.  (20)  is  small  compared  to  the 
linear  term.  In  this  case,  the  linear  and  nonlinear  terms  would  lie 
close  together. 

In  one -dimensional  blooming,  we  have  found  the  values  of 


aan(p) 

apk 


to  be  typically  in  the  range  from  0.  01  to  0.  05,  yielding  comparatively 
small  nonlinear  terms.  Thus,  one  hypothesis  on  the  superiority  of  the 
1-D  correction  is  that  the  polynomial  interactions  might  be  stronger  in 
a 2-D  than  in  a 1-D  medium. 

A second  hypothesis  is  that  the  nonlinear  terms  are  similarly 
r.mall  in  one  and  two  dimensions,  but  that  the  larger  number  of  degrees 
of  freedom  increases  the  nonlinear  correction  contribution.  It  we  were 


31 


to  expand  the  list  of  Legendre  and  Zernike  equivalents  (Section  2.  A.  2),  we 
would  find  that  if  a 1-D  blooming  medium  required  an  N-polynomial  descrip- 
tion, a 2-D  medium  with  comparable  spatial  frequencies  would  require 
approximately  (N^  + 4N)/4  terms. 

c . Predicting  the  Steady- State  Linear  Solution 

In  a linear  correction  scheme,  there  are  two  opposing  factors: 
linear  correction  drives  the  intensity  to  a peak,  but  medium  response  drives 
toward  a steady  state,  usually  diminishing  the  glint  intensity.  The  simula- 
tions we  ran  seemed  to  indicate  that,  within  limits,  the  relative  time  scales 
of  the  two  processes  are  unimportant.  This  encouraged  us  to  try  a model 
in  which  the  two  were  separated  completely:  a sequence  of  purely  linear 
correction  to  some  converged  value,  followed  by  the  determination  of  a new 
steady-state  medium.  When  the  linear  correction  holds  the  system  in  the 
same  state  as  the  previous  iteration,  we  can  assume  that  a steady- state 
system  solution  has  been  reached. 

Figure  5 shows  how  this  works  for  one-channel  tilt  correction.  The 
solid  line  gives  target  intensity  at  steady  state  as  the  transmitter  tilts  the 
focal  pattern  across  the  target.  The  maximum  of  this  curve  (at  |3j  = -19, 

IQ  = 0.43)  is  the  nonlinear  solution.  We  start  the  predictive  linear  model 
at  various  initial  tilt  positions.  The  arrows  indicate  the  direction  of  the 
system  in  converging  to  an  instantaneous  linear  correction.  The  dashed 
lines  give  the  change  in  intensity  as  the  medium  drives  toward  steady  state. 
These  sequences  converge  to  a range  of  potential  linear  solutions  about 
Pj  = 16,  IQ  = 0.39.  Only  a few  iterations  were  necessary  to  find  this  solu- 
tion, which  resulted  in  a huge  saving  of  computer  time. 

This  model  was  tested  on  the  same  cases  used  in  the  earlier  multi- 
dither simulations: 

• Case  I:  In  very  few  iterations,  the  system  converged  to  the 
same  point  as  the  earlier  simulations. 

• Case  II:  The  system  oscillated  initially  between  the  same 
limits  seen  in  the  full-scale  simulations,  then  settled  down 
to  a good  correction.  This  correction  was  then  used  as  an 
initial  condition  in  the  full  time -dependent  simulation.  The 
results,  shown  in  Figure  6,  indicate  that  the  predictive 
model  is  valid. 
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Figure  5.  Multidither  predictive  model  — one  channel. 
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Figure  6.  Multidither  predictive  model  results  — case  II. 
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The  nonlinear  and  linear  solutions  to  Case  II  are  presented  in 
Table  3.  We  see  that  they  are  again  very  similar,  as  predicted. 

Table  3.  Nonlinear  and  Linear  Results 


Steady-State 
Corrected 
Target  Intensity 


Nonlinear 

Linear  prediction 

Multidither 
starting  from 
linear  predic- 
tion results 


Polynomial  Coefficients 
(radians) 


-11  0.23  2.0  -1.8  -3.4 

-11  0.  32  2.0  -1.  1 -3.0 

-10.7  0.26  2.2  -1.3  -2.8 


Full-scale  simulations  in  2-D  have  been  too  expensive  to  permit  their 
steady- state  behavior  to  be  investigated  extensively.  However,  with  the 
predictive  linear  model,  this  becomes  feasible. 

As  indicated  above,  we  might  expect  to  see  a greater  discrepancy 
between  the  nonlinear  and  linear  steady-state  solutions  in  2-D.  We  checked 


this  fo 


a sample  case: 

Power  = 375  kW 
Wind  velocity  = 12.5  m/sec 
Focal  length  = 2 km 
a = 2 x 10"^/cm 
U = 0.  02  rdn/sec 


p =35  cm 
ro 

P,  , - 26.  75  cm  . 

l / e 


We  obtained  the  following  results: 
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Table  4.  Sample  Case  Results 


I 


Model 

On-Axis 

Intensity 

Zernike  Polynomial  Coefficients 
(radians) 

Z11 

Z20 

Z40 

Z31 

Z22 

Z33 

Z42 

Nonlinear 

0.42 

4.  0 

-0.  25 

-0.  12 

-1.  5 

1.6 

1.  2 

-0.  12 

Linear  prediction 

0.  33 

m 

-0.  98 

0.  01 

-1.4 

2.  8 

1.5 

-0.  87 

The  discrepancies  are  indeed  greater  and  the  correction  not  as  good  as 
in  1-D. 

C.  SEQUENCING  THE  POLYNOMIAL  CHANNELS 

1 . Multidither 

We  investigated  whether  changing  the  polynomial  sequences  would  create 
a maximization  path  from  an  initially  flat  mirror  to  a better  solution. 

Figure  7 shows  the  results  of  three  of  the  five  cases  which  were  run.  The 
arrows  indicate  the  time  at  which  a subsequent  channel  was  added,  and  the 
numbers  give  the  order  of  the  added  polynomial.  In  all  cases,  the  tilt  cor- 
rection was  the  first  to  be  stabilized. 

In  every  case,  the  sequencing  system  converged  to  the  same  secondary 
solution  found  by  the  simultaneous  channel  multidither  system.  In  this 
solution,  P5(x)  is  so  strongly  dominant  that  no  significant  correction  is  seen 
until  its  channel  is  cut  in. 

2.  Function  Maximization 

Our  FM  code  does  a form  of  sequencing  as  it  varies  one  polynomial 
coefficient  at  a time.  We  changed  the  sequence  of  the  polynomials  and  found 
that  the  path  could  be  shifted  from  one  maximum  to  another,  as  shown  in 
Table  5 for  case  I (the  FM  system). 


1.0 


6964-23 


Figure  7.  Multidither  Legendre  sequencing  — case  I. 
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Table  5.  Function  Maximization  Sequencing  Results 


Polynomial  Sequence 

Intensity 

Polynomial  Coefficients 

P2 

P4 

P5 

3 - 5 - 2 - 4 
2 - 3 - 4 - 5 

0.  895 

0.  095 

6.  97 

-0.  1 

-0.  5 

5 - 3 - 4 - 2 
4 - 2 - 5 - 3 
2 - 4 - 3 - 5 

0.  y7 

1 

o 

c 

0.  68 

-0.  22 

-3.  85 

We  do  not  know  whether  this  is  caused  by  the  nonlinearity  of  the  FM  solution 
or  by  the  freezing  of  previously  maximized  coefficients.  The  multidither 
sequencing  code  permitted  previous  channels  to  continue  optimising  as  new 
ones  were  added. 
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D. 


CONCLUSIONS 


We  found  that  correction  was  excellent  for  all  types  of  one -dimensional 
systems.  Function  maximization  (nonlinear  correction)  did  best,  followed 
by  PPC  (linear  correction)  and  multidither.  A comparison  of  temporal 
behavior  showed  that  fast  multidither  and  PPC  were  practically  equivalent. 
The  modal  1-D  multidither  solution  was  not  changed  by  servo  bandwidth 
variation  or  by  polynomial  sequence.  Also,  the  linear  and  nonlinear  solu- 
tions were  closer  together  in  1-D  than  in  2-D,  which  coupled  to  the  better 
correction  in  1-D  than  in  2-D. 


SECTION  3 


I 
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EXPERIMENTAL  BLOOMING  COMPENSATION 

This  section  presents  a detailed  discussion  of  the  experimental 
measurements  related  to  thermal  blooming  compensation  which  were  made 
during  this  contract.  Various  algorithms  are  described,  a 19- 
actuator  zonal  compensation  algorithm,  a modal  algorithm 
employing  Zernike  polynomials,  and  in  addition,  the  use  of  an 
auxiliary  focus  and  tracking  servo  is  discussed. 

A.  EXPERIMENTAL  SETUP 

The  experimental  measurements  were  made  on  the  DARPA/RADC 
COAT  system  modified  to  incorporate  the  deformable  mirror  and  the  nec- 
essary electronics  for  the  Zernike  polynomial  studies.  A detailed  descrip- 
tion of  these  modifications  is  given  in  interim  reports  2 and  3 and  will  not 
be  repeated  here.  ' 

The  optical  arrangement  used  in  making  these  measurements  is  shown 
in  Figure  8.  There  are  two  optical  paths,  one  known  as  the  "local  loop," 
which  is  an  undistorted  path,  and  a propagation  path  where  the  blooming 
distortions  are  introduced.  In  these  experiments,  a liquid  absorption  cell 
was  used  to  simulate  the  scaled  thermal  blooming  conditions  encountered 
by  focused  beams  propagating  over  relatively  short,  low-altitude  paths  in 
the  atmosphere.  This  cell  was  previously  used  on  contract  F30602-75-C- 
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0001  and  is  discussed  in  the  final  report  on  that  contract.  The  cell  (shown 
in  Figure  9)  is  20  cm  long  with  flow  transverse  to  the  length  dimension. 
Methanol,  with  a small  amount  of  iodine  added  to  absorb  the  0.488  pm  laser 
light,  was  the  liquid  used. 

In  a blooming  experimental  sequence  photomultiplier  (PMT)  1 is  con- 
nected to  the  COAT  servo  system  and  a near  diffraction- limited  beam  is 
produced  in  the  local  loop  at  target  plane  1.  In  addition,  this  beam  is  put 
into  the  propagation  path  containing  the  auxiliary  tracking  and  focus  controls 
and  the  liquid  blooming  cell.  The  peak  target  irradiance  (or  equivalently 
the  Strehl  ratio)  at  target  plane  2 is  then  measured.  This  gives  the  NO  COAT 


39 


4640  1 


PUMP 


. " 

OPTICAL 

BEAM 


GLASS  WINDOW 


FLOW 

DIFFUSING 

MANIFOLD 


TRANSVERSE  FLOW  CHANNEL 


METAL  AND  GLASS  BOX 


Figure  9.  Blooming  cell. 
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corrected  Strehl  ratio  of  the  bloomed  beam.  When  PMT  2 is  subsequently 
connected  to  the  COAT  servo  system,  the  degree  of  COAT  correction  of  the 
bloomed  beam  is  determined. 

The  diagnostics  used  to  analyze  system  performance  depended  on  the 
particular  experiment.  In  general,  it  consisted  of  TV  monitoring  and  pin- 
hole photodiode  detectors.  As  discussed  below,  the  bloomed  beam  shifts 
into  the  wind  and  the  peak  target  irradiance  does  not  in  general  occur  at  the 
boresight.  The  method  for  determining  the  peak  target  irradiance  for  the 
NO  COAT  corrected  bloomed  beam  consisted  of  translating  the  pinhole 
detector  to  the  point  in  the  beam  at  which  the  intensity  was  a maximum.  A 
second  method  which  was  also  used  consisted  basically  of  fixing  the  pinhole 
detector  and  moving  the  bloomed  beam  to  the  detector.  This  was  done  by 
using  the  auxiliary  tracking  mirrors  (see  below)  and  a servo  system  separate 
from  the  deformable  mirror.  The  two  methods  gave  the  same  results,  but 
the  later  method  was  operationally  simpler. 

F or  all  the  experiments  to  be  discussed  below,  the  following  termin- 
ology will  be  used,  except  where  noted:  "COAT  OFF"  means  PMT  1 is  used 
and  "COAT  CORRECTED"  means  PMT  2 is  used.  To  facilitate  comparing 
the  various  algorithms  studied,  the  local  loop  Strehl  ratio  was  monitored 
continuously  throughout  the  experiments  with  the  pinhole  detector  shown  in 
the  local  loop  path  in  Figure  8. 

B.  ALGORITHMS 

Within  the  limitations  of  the  experimental  equipment,  various  temporal 
and  spatial  algorithms  were  investigated  to  determine  their  degree  of  com- 
pensation for  thermal  blooming.  Included  for  spatial  algorithms  were  a 
19-actuator  deformable  mirror  zonal  mode  and  a 5-Zernike  polynomial 
modal  mode.  For  temporal  algorithms,  sequencing  of  the  5 Zernikes  was 
included.  In  addition,  the  auxiliary  tracking  and  focus  controls  were  studied 
with  various  combinations  of  the  above  algorithms.  The  tracking  and  focus 
servos  have  dither  frequencies  of  1 kHz  and  100  Hz,  respectively.  It  was 
hoped  that  these  slower  dithers  would  represent  a temporal  algorithm 
different  from  the  deformable  mirror  system  at  the  higher  dither  frequencies 
(~10  to  30  kHz).  Consequently,  track  and  focus  represent  both  a temporal 


42 


and  a spatial  algorithm.  From  the  analytic  results  (see  Section  2),  it 
appears  that  dither  frequencies  much  lower  than  100  Hz  are  required  before 
a temporal  effect  on  blooming  compensation  can  be  observed.  This  is  con- 
sistent with  the  observed  experimental  results  (see  below). 

C.  SCALING  CONSIDERATIONS 

The  laboratory  experiments  performed  on  this  contract  are  scaled  so 
that  the  results  can  be  extended  to  other  wavelengths,  ranges,  etc.  In 
particular,  our  visible  wavelength  experiments  are  chosen  to  scale  to  a 
typical  10.6  pm  scenario. 

5 

According  to  Bradley  and  Herrmann,  the  propagation  of  a slewed 
beam  through  an  absorbing  medium  can  be  characterized  by  four  dimension- 
less parameters: 


Absorption  number  N^ 


= oR 


Fresnel  number  Np  = Ka  /R 


Slewing  number  N 


Distortion  number  N, 


4Tm,j,aRP 
npCpVaX.  ’ 


where  a is  the  absorption  coefficient,  R is  the  range,  K is  the  wave- 
number,  a is  the  l/e  radius  at  the  cell  entrance,  u is  the  slewing  frequency, 
v is  the  cross-wind  velocity,  n^,  is  the  coefficient  of  index  change  with 
respect  to  temperature,  n is  the  index  of  refraction,  P is  incident  laser 
power,  p is  medium  density,  Cp  is  specific  heat  at  constant  pressure  of  the 
medium,  and  X is  wavelength  of  the  laser.  In  addition,  for  typical  10.6  pm 
thermal  blooming  scenarios,  heat  transfer  in  the  beam  is  convection  rather 
than  conduction  dominated.  Hence,  for  laboratory  experiments,  the  wind 
speed  and  beam  size  must  be  chosen  to  ensure  a convection-dominated 
situation.  This  can  be  satisfied  by  requiring  that  the  Peclet  number  be  much 
greater  than  one.  The  Peclet  number,  which  is  the  ratio  of  convection  to 
conduction,  is  given  by 

„ 2va 
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where  K is  the  medium  thermal  conductivity.  For  the  present  experiments, 
Pe  ~ 50;  hence,  convection  dominates,  the  distortion  number  varies  from 
0.6  to  60,  the  Fresnel  number  is  typically  130,  the  absorption  number  is 
~0.2,  and  there  is  no  slewing. 

Although  these  parameters  characterize  the  blooming  medium,  it 
is  crucial  to  know  the  relationship  of  the  focused  beam  path  to  the  length  of 
the  blooming  cell  and  to  the  position  of  the  blooming  cell  in  the  beam  path. 
In  previous  studies,  the  effect  of  near-field  or  far-field  blooming  was 
presented,  illustrating  this  importance.  In  the  present  laboratory  experi- 
ments, the  entire  focused  beam  path  is  within  the  blooming  medium  as  it 
would  be  in  a scaled  scenario. 


D.  ZONAL  COAT  PERFORMANCE 

A typical  intensity  profile  for  both  a local  loop  beam  and  a bloomed 
beam  with  a high  distortion  number  (ND  ~ 20)  is  shown  in  Figure  10. 
Qualitative  comparisons  in  amplitude  should  only  be  inferred  because  the  TV 
camera  used  to  generate  this  data  is  nonlinear.  This  data  is  actually  taken 
with  the  Zernike  mode,  but  in  fact  the  data  in  the  19  actuator  zonal  mode  is 
essentially  the  same. 

Figure  11  gives  a sequence  of  COAT  OFF  and  COAT  CORRECTED 
photographs  of  the  beam  profile  taken  for  three  distortion  numbers,  from  a 
condition  of  essentially  no  blooming  to  a case  of  strong  blooming  (the  f#  of 
the  TV  camera  is  shown).  We  see  the  expected  crescent- shaped  beam  for 
the  strong  blooming  case.  Further,  the  beam  was  shifted  into  the  wind,  as 
expected  (not  shown).  Figure  12  shows  the  peak  target  irradiance  as  a func- 
tion of  input  power  for  the  data  of  Figure  11.  These  results  are  for  N.  =0.2. 

A 

Np  = 318  P,  and  Np  = 134.  Plotted  are  results  for  COAT  OFF,  COAT  COR- 
RECTED, and  for  no  blooming  (linear  with  power).  The  COAT  system 
increases  the  peak  target  irradiance  by  up  to  a factor  of  1 . 8 at  a power  input 
of  ~20  mW.  But  as  the  power  is  further  increased,  correction  decreases 
until  virtually  no  correction  is  achieved  at  the  higher  powers.  However,  we 
can  see  from  Figure  11  that,  at  the  high  power  condition,  there  is  a correc- 
tion obtained  which,  although  it  would  not  affect  the  peak  target  irradiance, 
would  result  in  a higher  irradiance  if  a larger  aperture  were  used. 


44 


5964-1 1 


INTENSITY 

(BEAM  INTENSITY  5 x LESS 
THAN  BLOOMED  PATH  BEAM) 


LOCAL  LOOP 
COAT -ON 


INTENSITY 


BLOOMED  PATH 
NO  COAT  CORRECTION 
P = 56  mW 


Figure  10.  Intensity  profiles. 
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Figure  11.  Photographs  of  bloomed  beam 
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Absolute  Strehl  ratio  was  determined  for  this  case  and  the  results  are 
shown  in  Figure  13.  In  determining  the  Strehl  ratio,  a truncated  gaussian 
near-field  input  beam  profile  of  measured  truncation  radius  and  l/e  inten- 
sity was  used.  The  radius  measurements  were  made  using  the  scanning 
technique  discussed  in  interim  report  3.  The  no-blooming  case  Strehl  ratio 
is  about  80%.  We  see  that,  when  the  bloomed  Strehl  ratio  with  COAT  OFF 
drops  to  about  30%,  the  COAT  system's  ability  to  continue  producing  1 .8 
correction  factors  ceases. 

The  above  results  were  achieved  with  the  19-actuator  zonal  mirror. 

The  auxiliary  tracking  and  focus  controls  were  also  used.  When  using  the 
tracking  control  and  the  deformable  mirror,  the  peak  target  irradiance  was 
not  changed.  This  is  expected  based  on  the  method  of  measuring  the  peak 
target  irradiance  described  previously.  That  is,  the  pinhole  detector  is 
translated  to  the  point  of  maximum  intensity  on  the  beam.  Since  the  tracking 
control  acts  the  same  as  translating  the  pinhole  detector,  no  difference  is 
observed.  However,  fcr  determining  boresight  irradiance,  tracking  control 
is  essential  and,  in  this  case,  keeps  the  beam  on  the  boresight  (or  to  wher- 
ever the  target  glint  is  moved).  It  was  not  possible  to  use  the  auxiliary  focus 
control  alone  because  of  spherical  distortion  in  the  focus  actuator  (see 
Ref.  4 for  a description  of  this  problem).  However,  the  auxiliary  focus 
control  was  used  in  conjunction  with  the  zonal  mirror.  In  this  case,  no 

A 

additional  correction  was  observed.  This  is  consistent  with  previous  work 
and  is  expected  based  on  the  dither  results  discussed  in  the  analytic  section 
(Section  2). 

Figures  14,  15,  and  16  give  similar  plots  of  peak  target  irradiance  as 
a function  of  input  power  for  three  sets  of  the  parameters  N , N_,  and  Is'  . 

A U r 

Correction  was  achieved  in  all  cases.  Of  the  three  parameters,  it  appears 
that  the  Fresnel  number  is  the  most  critical  in  establishing  the  level  of 
blooming  and  the  level  of  COAT  correction,  with  the  best  correction  achieved 
for  the  cases  with  large  Fresnel  numbers.  For  these  cases,  operation  with 
the  auxiliary  focus  and  tracking  controls  was  also  performed  with  results 
consistent  with  those  discussed  above. 

One  test  was  performed  with  input  power  fixed  but  wind  speed  variable. 
The  results  of  this  test  are  shown  in  Figure  17,  which  gives  peak  target 
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Figure  14.  Peak  target  irradiance  versus  input  power  — 
zonal.  Na  = 0.2,  Np  = 200P,  Np  = 340. 
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Figure  16.  Peak  target  irradiance  versus  input  power  — 
zonal.  Na  = 0.8,  ND  = 1270P,  Np  = 134. 
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irradiance  as  a function  of  wind  speed.  Intensity  increases  nearly  linearly 
over  the  range  of  the  experimental  conditions.  This  is  the  expected  trend 
based  on  the  argument  that  distortion  number  is  inversely  proportional  to 
wind  speed. 

E.  MODAL  COAT  PERFORMANCE 

As  described  in  Ref.  1 and  2,  seven  Zernike  polynomials  were  wired 
into  the  DARPA/RADC  COAT  system.  These  are  shown  in  Table  6. 


Table  6.  Zernikes  in  DARPA/RADC  System 


Zernike 

Function 

Error  (%) 

1 

Horizontal  tilt 

7 

2 

Vertical  tilt 

7 

3 

Astigmatism 

13 

4 

Astigmatism 

13 

5 

F ocus 

23 

6 

Coma  (y- symmetry) 

51 

7 

Coma  (x- symmetry) 

51 

The  deformable  mirror  was  made  to  conform  approximately  to  these  Zernike 
functions  through  the  use  of  a set  of  numbers  representing  displacements  of 
the  various  actuators.  This  set  was  generated  by  a function-maximization 
routine,  the  details  of  which  are  discussed  in  Section  4.  As  shown  in 
Section  4,  large  errors  exist  between  the  actual  analytic  Zernike  function 
and  the  approximation  to  it  made  by  a finite  actuator  mirror,  particularly 
for  the  higher-order  functions.  For  the  19-actuator  mirror  used  in  the 
present  experiments,  the  errors  are  given  above.  Large  errors  in  the 
functions  imply  that  they  are  not  orthogonal  and,  consequently,  that  cross 
coupling  would  be  expected  to  affect  the  results.  Modal  operation  of  the 
COAT  system  in  the  unbloomed  (local-loop)  path  gave  evidence  of  this.  To 
investigate  servo  stability,  each  of  the  seven  Zernikes  was  switched  indi- 
vidually into  the  COAT  system.  Each  one  converged  to  a stable  state. 
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However,  when  all  were  switched  on  sequentially,  an  unstable  state  resulted 
when  the  higher  order  Zernikes,  Z6  and  Z7,  were  included.  That  is,  with 
Z1  through  Z5  on,  switching  Z6  and/or  Z7  on  resulted  in  an  unstable  servo. 
This  probably  resulted  from  the  large  50%  error  in  the  functions  and  the 
cross-coupling  caused  by  the  error.  In  subsequent  experiments  concerning 
blooming  compensation,  the  results  from  which  are  discussed  below,  only 
the  first  5 Zernikes  listed  were  used. 

Experiments  performed  for  blooming  compensation  used  three  sets  of 
the  parameters  N^,  N^,  and  Np.  The  results  of  these  experiments  are 
shown  in  Figures  18,  19,  and  20.  These  show  that  in  only  one  case  (for 
= 0*2,  Kp  = 318  P,  and  Np  = 134,  as  shown  in  Figure  18)  did  modal 
COAT  using  the  5 Zernikes  give  a correction.  Comparing  this  result  to  the 
corresponding  zonal  case  (Figure  12)  indicates  a reasonably  similar  result 
although  the  correction  is  somewhat  poorer.  However,  the  Zernike  mode 
gives  no  apparent  correction  for  the  other  cases.  Again,  using  the  auxiliary 
focus  and  tracking  servos  yielded  the  same  result  as  had  the  zonal  case. 

It  is  important,  when  comparing  the  zonal  and  modal  cases,  that  the 
local  loop  beam  profile  be  the  same  in  each  case.  This  is  so  because  it  is 
this  beam  that  is  propagated  into  the  blooming  cell.  Strehl-ratio  measure- 
ments of  the  local  loop  beam  in  the  modal  mode  gave  virtually  the  same 
result  as  those  in  the  zonal  mode.  Comparing  the  results  of  Figure  12  with 
those  of  Figure  18  shows  that  the  COAT  OFF  results  are  slightly  better  in 
the  modal  mode.  Comparing  Figure  17  with  Figure  19  and  Figure  16  with 
Figure  20  shows  that  the  modal  COAT  OFF  result  is  nearly  the  same  as  the 
zonal  corrected  result.  This  implies  that  the  near-field  beam  input  to  the 
cell  may  be  different.  One  would  expect  that  the  corrections  would  be  better 
if  the  blooming  were  not  as  severe.  Unfortunately,  this  was  not  borne  out 
by  the  results.  This  may  have  been  because  of  the  nonorthogonality  of  the 
Zernikes  or  because  higher  order  Zernikes  would  have  been  needed  to  give  a 
correction.  Earlier  analytic  work  on  the  contract  showed,  basically,  that 
seven  select  Zernikes  were  necessary  to  give  good  correction  in  a moderate 
thermal  blooming  case.  These  are  horizontal  tilt,  focus,  astigmatism 
(x  - y only),  coma  (y-symmetry  only),  spherical  aberration,  5th  order 
coma,  and  5th  order  spherical  aberration.  Thus,  for  the  present 
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experimental  measurements,  only  horizontal  tilt,  focus,  and  coma  are 
effective.  The  remaining  four  functions  are  higher  order  and  would  require 
more  actuators  than  the  present  19. 

Returning  to  the  case  where  correction  was  achieved,  a sequencing  set 
of  experiments  with  the  five  Zernikes  was  performed.  The  tabulated  data 
for  one  power  level  for  various  combinations  of  Zernikes  is  shown  in 
Table  7.  Although  other  cases  were  run,  these  results  are  typical.  Power 
was  0.42  (run  1)  with  all  five  on.  Power  dropped  to  0.32  when  Z5  was 
turned  off  for  run  2.  Power  dropped  to  0.34  when  Z4  was  turned  off  for 
run  3.  No  power  change  was  observed  for  runs  4,  5,  and  6.  Runs  4,  5, 
and  6 correspond  to  astigmatism  (x,*  y),  y tilt,  and  x tilt.  With  all  off 
(run  7),  power  was  0.24.  These  results  are  consistent  with  the  analytic 
results  from  Ref.  4 discussed  above  and  with  the  computer  simulation  work 

Table  7.  Zernike  Sequencing 


Sequence 


all-on 


Peak  Irradiance  (arb.units) 


all  off 


54321 


34512 


other 


2 2 

discussed  in  Section  2.  That  is,  focus  and  astigmatism  (x  - y ) are  the 
most  important  functions  in  compensating  for  blooming.  The  fact  that  tilt 
was  not  found  to  be  important  was  a result  of  the  measurement  technique 
(described  previously).  That  is,  since  the  pinhole  detector  was  always 
translated  to  the  point  of  maximum  intensity,  tilt  could  not  affect  measured 
peak  intensity. 

Runs  8,  9,  and  10  represent  sequences  in  which  the  5 Zernikes  were 
turned  on.  No  difference  in  peak  irradiance  was  obtained  for  any  sequence 
(run  10  represents  a large  number  of  other  sequences).  What  we  did  observe 
was  that,  regardless  of  when  Z4  and  Z5  were  turned  on  in  the  sequence,  the 
final  result  was  always  the  same. 


60 


SECTION  4 


I 

4 


zernike-polynomial/deformable-mirror 

APPROXIMATION 

The  Zernike  polynomials  and  the  calculation  procedure  used  to  generate 
the  actuator  inputs  (which  were  required  to  make  the  deformable  mirror 
approximate  those  polynomials)  are  described  in  this  section.  This  set  of 
actuator  inputs  was  implemented  in  the  DARPA/RADC  deformable  mirror, 
as  discussed  in  Section  3. 

A.  ZERNIKE  POLYNOMIALS 

The  Zernike  polynomials  are  a complete  set  of  orthogonal  functions 
defined  on  the  unit  circle.  A useful  table  of  the  Zernike  polynomials  in  polar 
form  is  given  in  the  Perkin-Elmer  report  RADC-TR-76- 1 16.  The  first  14 
of  these  functions  and  their  equivalents  in  rectantular  coordinates  are  repro- 
duced in  Table  8.  The  seven  implemented  in  the  DARPA/RADC  electronics 
(as  described  in  Section  3)  were  Z1  through  Z7. 

B.  PROBLEM  DEFINITION 

A deformable  mirror  is  typically  a thin,  flat,  optically  polished  plate 
supported  on  the  back  by  an  array  of  extensible  actuators.  The  actuators 
are  attached  to  a massive,  rigid  baseplate.  The  actuators  are  usually 
cylinders  or  disc  stacks  of  piezoelectric  ceramic  which  can  be  controlled 
electrically  to  produce  the  desired  extensions  and  thus  to  deform  the  mirror 
plate  in  a controlled  manner.  The  total  travel  of  an  actuator  is  limited  to  a 
few  wavelengths  of  the  laser  radiation  so  the  deformation  of  the  plate  remains 
within  the  elastic  limit,  and  the  linear  theory  of  plates  is  applicable. 
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Ze  rnike 

Polar  Constant 

Rectangular 

1 

2 r cos  6 1 

2 x 

l Tilt  (lateral) 

2 

2 r sin  6 j 

* y 

3 

\rb  r2  sin  26  1 

1 astigmatism 

2 6 xy 

4 

sTb  r2  cos  26  J (3rd  °rder) 

\"6  (x2  - y2) 

c 

\3  (2 r2  - 1)  defocus 

sri  (2(x2  + y2)  - 1) 

(longitudinal) 

6 

\f8  (3r3  - 2r)  sin  6 1 

1 coma 

\r8  y (3(x2  *■  y2)  - 2) 

pi 

vr8(3r3  - 2r)  cos  6 J *3rd  °rder> 

sIS  x (3(x2  + y2)  2) 

IS 

nT8  r3  .in  36  1 ,20o 

sTS  y (3x“  - y2) 

H 

r3  co.  36  J (fifth  order> 

sf B x (x2  - 3y2) 

M 

■v/5  (6r4  - 6r2  * 1)  3rd  order  spher. 

sfe  (6(x2  + y2)2  - 6(x2  + y2)  1) 

m 

vlon/.s,2,  co.  1.1 

sTlO  (x2  - y2)  (4  (x2  + y2)  - 3) 

m 

«^0  (4r4  - 3r2)  .in  26  J 5th  ord*r 

2\r10  xy(4(x2  + y2)  - 3) 

m 

sTlO  r4  co.  46 

sflO  ((x2  * y2)2  - 8x2y2) 

s/lO  r4  .in  46 

4««/l0  xy(x2  - y2) 
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In  a linear  system,  the  principle  of  superposition  states  that 
the  plate  deformation  produced  by  all  the  actuators  operating  together 
can  be  expressed  as  the  algebraic  summation  of  the  deformations 
produced  by  the  individual  actuators  acting  separately.  For  each 
actuator,  a unit  profile  can  be  defined.  It  is  the  profile  produced  on 
the  mirror  when  that  actuator  receives  unit  electrical  input  and  all 
other  actuators  receive  zero  input.  The  total  profile  can  then  be 
calculated  as  the  point-by-point  summation  over  the  surface  of  all  » 

the  unit  profiles,  each  multiplied  by  the  algebraic  value  of  its  electrical 
input  at  that  instant. 

In  a real  mirror,  the  various  unit  profiles  will  not  be  identical 
because  of  distortion  produced  by  proximity  to  the  edge  of  the  mirror. 

In  this  study, these  distortions  will  be  neglected  and  it  is  assumed  that 
all  actuators  will  be  characterized  by  identical  unit  profiles. 

Laboratory  measurements  made  on  several  mirrors 5 show  that 
the  unit  profiles  can  be  characterized  with  good  accruacy  by  a modified 
"gaussian"  having  circular  symmetry  and  defined  as: 

Z = exp  (A  rn) 

The  true  gaussian  corresponds  to  a value  of  n = 2,  while  the  measured 
profiles  ranged  from  n = 1.5  to  n = 2.5.  The  sign  of  A is  negative  and 
its  numerical  value  determines  the  "width"  of  the  unit  profile  and 
hence  the  degree  of  cross  coupling  between  adjacent  actuators.  This 
cross  coupling  is  defined  as  the  motion  produced  at  the  center  of  gt 
given  actuator  when  an  adjacent  actuator  is  driven  with  unit  input. 

The  measured  values  have  ranged  from  3 to  20%,  In  particular,  for 
the  DARPA/RADC  mirror,  A = -5.369  and  n = 1.71  and  the  measured 
cross  coupling  was  14%. 
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C.  SELECTED  PARAMETER  VALUES 

For  this  study,  the  actuator  array  consisted  of  37  actuators 
arranged  in  a "circular-hex"  pattern  (as  shown  in  Figure  21).  In 
addition,  a 19-actuator  study  was  made.  This  was  done  by  not  using 
the  outer  ring  of  actuators  shown  in  Figure  21.  The  results  from  this 
19-actuator  study  were  implemented  into  the  DARPA/RADC  mirror. 
However,  to  facilitate  an  understanding  of  the  computational  scheme, 
the  37-actuator  system  will  be  described  in  detail  with  the  19  actuator 
results  summarized  briefly. 

The  unit  circle  of  the  Zernike  functions  passes  through  the 
centers  of  the  outer  circle  of  actuators.  Also  shown  in  Figure  21  is 
the  array  of  mesh  points  on  which  the  numerical  calculations  are  per- 
formed. The  mesh  points  also  include  a point  at  the  center  of  each 
actuator.  The  total  number  of  mesh  points  is  127.  The  expression 
for  the  unit  profile  for  the  37-actuator  study,  which  was  chosen  as  a 
true  gaussian  having  a coupling  of  12%  to  adjacent  actuators,  is 

z = exp  (-9.  0824  r2)  , 

where  r is  measured  from  the  center  of  the  given  actuator  and  the 
scale  was  chosen  so  that  the  unit  circle  has  a radius  of  1. 

D.  COMPUTER  PROGRAM 

The  goal  of  the  calculation  was  to  find  a vector  of  actuator 
inputs  such  that  the  form  of  the  mirror  would  be  optimally  matched 
to  the  form  of  a given  Zernike  polynomial  in  the  least-squares  sense. 
Several  methods  were  tried  which  failed  because  of  poor  accruacy  or 
excessive  CPU  time.  These  included  multivariate  optimization  (hill 
climbing)  techniques,  simulation  of  a complete  COAT  system  matching 
a phase  front  in  the  form  of  a Zernike  polynomial,  and  a method  of 
local  minimization  which  reduced  the  error  in  the  neighborhood  of 
each  actuator.  The  last  method  gave  good  results  but  was  not  a true 
least- squares  solution. 
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The  method  which  was  finally  adopted  involved  application  of 
a computer  library  routine  called  LLSQAR.  The  main  program  starts 
by  generating  a table  of  actuator  locations  and  a table  of  mesh  point 
locations,  both  in  x,  y coordinates.  Next,  the  desired  Zernike 
polynomial  is  defined  as  a table  of  displacements  at  each  of  the  127 
mesh  point  locations.  Then  a matrix  table  of  influence  coefficients 
is  constructed.  This  is  a 37  x 127  table  which  gives  the  value  of  the 
unit  profile  for  each  actuator  at  each  mesh  point  location.  These  tables 
are  then  used  to  define  a set  of  127  linear  equations  in  37  unknowns. 

The  individual  actuator  inputs  are  the  unknowns,  the  influence  coef- 
ficients comprise  the  coefficient  matrix  of  the  equations,  and  the  table 
of  desired  Zernike  deflections  at  the  mesh  points  becomes  the  right- 
hand  vector  of  the  equations.  These  equations  constitute  an  over- 
determined set  with  more  equations  than  unknowns. 

The  function  of  the  subroutine  LLSQAR  is  to  solve  an  over- 
determined set  of  linear  equations.  It  yields  a set  of  values  for  the 
unknowns  which  will  best  satisfy  the  right-hand  values  with  a minimum 
total  squared  error.  The  result,  then,  is  a set  of  actuator  inputs 
which  produce  a mirror  surface  which  is  a least- squares  best  fit  to 
the  given  Zernike  polynomial  at  the  127  mesh  points. 

E.  OUTPUT 

The  results  of  the  calculation  are  given  as  tables  which  show 
the  actuator  inputs  required  to  approximate  each  of  the  Zernike  poly- 
nomials. In  addition,  a perspective  plotting  routine  was  used  to  pro- 
duce 3-D  displays  of  each  Zernike  function,  of  the  corresponding 
mirror  profile,  and  of  the  error  (which  is  the  difference  between  the 
profile  and  the  Zernike  function).  For  each  of  these  displays,  the 
data  is  recomputed  on  a 41  x 41  square  matrix  with  all  points  outside 
the  unit  circle  being  set  to  zero.  Finally,  the  rms  error  of  the  approxi- 
mation is  computed  and  printed.  This  number  is  derived  from  the 
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41  x 41  error  matrix  and  is  computed  as  the  square  root  of  the  mean 
squared  value  in  this  matrix,  again  disregarding  those  points  outside  the 
unit  circle.  The  matrix  contains  41  x 41  =1681  points  of  which  1245 
lie  within  the  unit  circle. 

Each  Zernike  polynomial  is  normalized  so  that  its  mean  square  value 
is  unity,  which  means  that  the  computed  errors  can  be  read  directly  as  per 
unit  values.  These  values  are  shown  on  the  corresponding  displays. 

A typical  result  for  the  actuator  inputs  for  the  37 -actuator  case  for 
the  focus  Zernike  is  given  in  Table  9.  Representative  perspective  plots 
for  the  tilt,  focus,  and  coma  Zernike  functions  are  shown  in  Figures  22, 

23,  and  24,  respectively.  We  see  that,  in  general,  generating  higher  order 
functions  with  the  mirrors  produces  larger  errors. 

In  Table  10,  the  rms  errors  for  the  37  actuators  and  the  19-actuator 
mirror  are  given  for  all  14  Zernike  functions.  It  is  apparent  from  the 
errors  shown  in  this  table  that,  to  generate  higher  order  functions 
accurately,  requires  more  actuators. 
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Table  9.  Focus  Zernike  — 37  Actuato 


1.  -O.  9454534 

2,  -0.94  54 *53  I 
.3,  -0.94  54  532 

4,  -0.9454530 

5,  -0.9454531 

6,  -0.9454530 

7,  -0.3925271 
5,  -0.392  52  70 

9,  -0.392  5270 

10,  -O'.  392  5271 

11,  -0.392  5272 

12,  -O'.  392  52  70 

13,  -O'.  39252  74 

14,  -O'.  39252  73 

15,  -0.392527  1 
15,  -0.3925270 
17,  -0.3925274 
J8,  -0.3925274 

1 9 , 1.4  52  0 72 

20,  I '.4  53000 

21,  r.  4 52072 

22,  1.4  52072 

23,  1'.  4 53508 

24,  1.4-52072 

25,  1. 4 52072 

26,  I'.  453898 

27,  I'.  452072 

28,  1.452072 

29,  1 .453808 

30,  1 .452072 

31,  1.452072 

32,  1.4  538  00 

33,  1.4  52072 

34,  I'.  4 52072 

35,  I'.  4 53898 

36,  1.452  0 72 

37,  -I'.  072  329 

RMIS  ERROR  = 0. 
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Table  10.  Zernike  Polynomial  RMS  Errors 


Zernike 

Error  for  37  Actuators, 

% 

Error  for  19  Actuators, 

% 

1 

5.2 

6.6 

2 

5.2 

6.  6 

3 

8.5 

12.8 

4 

8.5 

12.8 

5 

13.  1 

22.7 

6 

26.7 

50.7 

7 

26.7 

50.7 

8 

12.7 

21.9 

9 

12.7 

21.9 

10 

55.  3 

95.2 

11 

43.3 

78.  8 

12 

43.3 

78.8 

13 

18.4 

- 

14 

18.4 

- 
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SECTION  5 


RECOMMENDATIONS  FOR  FUTURE  WORK 

Although  the  experimental  and  analytic  results  obtained  from  this  work 
indicate  definite  trends  as  to  what  algorithms  can  be  used  for  blooming 
compensation,  the  results  have  also  indicated  several  areas  for  further 
investigation.  In  particular,  to  achieve  a more  conclusive  comparison 
between  the  zonal  and  the  modal  systems  using  Zernike  polynomials  will 
require  more  actuators  in  the  deformable  mirror.  This  would  permit 
higher-order  Zernikes  to  be  used  and  prevent  a large  degree  of  cross 
coupling  for  those  used.  In  addition,  it  would  be  desirable  to  have  a mirror 
with  a more  consistent  surface  excursion  from  actuator  to  actuator  (see  the 
third  interim  report).  This  would  require  time-consuming  polishing  of  the 
PZT  stacks  for  each  actuator. 

The  additional  capabilities  of  slewing  the  beam  would  be  a natural 
extension  of  the  blooming  scenarios  studied  under  the  contract.  For  the 
numerical  studies,  the  most  important  extension  is  to  perform  the  same 
type  of  study  using  the  more  costly  2-D  routines  just  made  operational  at 
Hughes.  With  such  results,  a detailed  quantitative  comparison  of  theory 
and  experiment  could  be  made. 
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METRIC  SYST1M 


BASE  UNITS: 

ja»e»J!!j/ 


length 

IMM 

lima 

electric  currant 
thermodynamic  temperature 
amount  of  substance 
luminous  intansity 

SUPPLEMENTARY  UNITS: 
plane  angle 
•olid  angle 

DERIVED  UNITS: 
Acceleration 

activity  (of  a radioactive  source) 

angular  acceleration 

angular  velocity 

area 

density 

electric  capacitance 

electrical  conductance 

electric  field  strength 

electric  inductance 

electric  potential  difference 

electric  reei stance 

electromotive  force 

energy 

entropy 

force 

frequency 

Illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  of  electricity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velocity 

viecosity.  dynamic 

viscosity  kinematic 

voltage 

volume 

wavenumber 

work 

SI  PREFIXES: 


. Unit 

metre 

kilogram 

Ncond 

imptrt 

kalvin 

mole 

candela 


radian 

ateradian 


metre  per  second  squared 
disintegration  per  second 
radian  par  second  squared 
radian  per  second 
square  metre 
kilogram  per  cubic  metre 
farad 
siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  par  kelvin 

newton 

harts 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

joule 

watt  per  ateradian 
joule  per  kilogram-kelvin 

pascal 

watt  per  metre-kelvin 
metre  per  second 
pescal-sscood 
square  metre  per  second 
volt 

cubic  metre 
reciprocal  metre 
joule 


..Multiplication  Factors 

I 000  000  000  000  - 10" 

1 000  000000  « 10* 
1000000-  to* 
1000-  10’ 
100  - 10> 
10-  10' 
0.1  - 10-' 
0.01  - 10-» 
0.001  - 10-> 
0.000  001  - to-* 
0.000  000  001  - to-* 
0.000  000  000  001  - 10“  " 
0000  000  000  000  001  - 10-" 
0000  000  000  000  000  001  • 10-'* 


SI  Symbol  Formula 


kg 

s 

A 

K 

mol 

cd 


tad 


m/s 

(disintegration  )/s 

red/s 

red/s 


F 

S 

H 

V 

V 

I 

N 

Hi 

lx 

Im 

Wb 

T 

A 

W 

Pa 

C 

I 


Pa 


V 

i 


kg/m 

AW 

A/V 

V/m 

Vs/A 

W/A 

VIA 

W/A 

N-m 

J* 

kg- m/s 
jcycleVs 

cd/m 
Oder 
Aim 
V-s 
Wb Im 

i/» 

Him 


W/sr 

j/kg-K 


W/m-K 

m/s 

Pee 


m/s 

W/A 


m 

(wavej/m 

N-m 


Prefix  81  Symbol 

tare 
SIB* 

kilT 

hecto* 
delu- 
ded* 
cent  I* 
mill! 
micro 


•NO 


To  he  avoided  where  possible 


nr*  aovtseusin  raarrmo  ornot:  — tii  «si  is 
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